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ABSTRACT 

We analyze the feasibility to estimate the stellar mass of galaxies by mid- 
infrared luminosities based on a large sample of galaxies cross-identified from 
Spitzer SWIRE fields and SDSS spectrographic survey. We derived the formu- 
lae to calculate the stellar mass by using IRAC 3.6/im and 4.5//m luminosities. 
The mass-to-luminosity ratios of IRAC 3.6/im and 4.5/rni luminosities are more 
sensitive to star formation history of galaxies than other factors, such as the in- 
trinsic extinction, metallicity and star formation rate. To remove the affection 
from star formation history, we used g-r color to recalibrate the formulae and ob- 
tain a better result. It must be more careful to estimate the stellar mass of low 
metallicity galaxies using our formulae. Due to the emission from dust heated 
by hottest young stars, luminous infrared galaxies present higher IRAC 4.5 /im 
luminosity compared to IRAC 3.6 /im luminosity. For most of type-II AGNs, the 
nuclear activity can not enhance 3.6/im and 4.5/im luminosities compared with 
normal galaxies. The star formation in our AGN-hosting galaxies is also very 
weak, almost all of which are early-type galaxies. 
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Introduction 



It is important to obtain the mass of galaxies for understanding the evolution of the 
universe. We know that the dark matter halo contains most of mass, and can only be 
detected using the effect of gravitation. So generally, the total mass of galaxies could only 
be computed on the basis of kinematics of bright stars, clusters, or even satellite galaxies (e.g., 
Rogstad fc Shostak Ill972j ; iRoberts fc Rots Ill973l ; lOstriker et al.lll974j ). However, due to the 
restriction from lower resolution and lower sensibility of telescopes for detecting the motion 
of celestial bodies in high redshift galaxies, these methods are just used for deriving general 
properties of nearby galaxies. Though the fraction of baryonic components in galaxies is 
relatively small, it could be detected and studied easily. Except the gas dominant galaxies, 
stars hold most of baryonic components in galaxies. Thus, to understand the correlation 
between the stellar mass and total mass is essential for studying the faint and distant galaxies. 
Therefore, it is necessary to seek some methods to estimate the stellar mass in galaxies. 

If we know the mass-to-luminosity ratio of a galaxy, and assume an initial mass function 
(IMF), its s tellar mass could also be estimated from the luminosity in corresponding wave- 
length (e.g.. lBell fc de Jong Il200ll ; iPortinari et al. II2004I ). Traditionally, mass-to- luminosity 
ratios could be derived by fitting photometric colors or spectra with models. Most of bary- 
onic components are trapped in low mass stars, which dominate the output of galaxies in 
longer wavelength range compared to the massive ones. At the same time, the short life 
scale of massive stars and heavy extinction in relatively short wavelength ranges restrict 
us to use them to detect the bulk of galactic stellar masses. Hence, less affected by other 
factors, such as metallicity, star formation rate (SFR) and star formation history (SFH), lu- 
minosities in relatively longer w avelength (for example, near-infrared (N IR) band), are ideal 



tracers of the stellar mass (e.g., ICole et al.ll2001t iKochanek et al.ll200ll ). The K band lumi- 



no sity distribution of dist a nt ga laxies was discussed in support of the hierarchical picture 
by iKauffmann fc Chariot I (119981 ) . The main uncertainty in the inferred stellar masses arises 



from the age of the stel lar population 



Rix fc Rieke 1119931 ). similar to the results calculated 



by optical luminosities (IBell et al.ll2005f). The discrepan cy among stellar masses derived from 



various methods has been found (e.g.. lDrory et al.ll2004l ). For example, th e uncertainty of the 



stellar mass estimated by if-band lum inosity alone has been discussed by iBrinchmann et al 



( 120001 ) . Kannappan &: Gawiser I (120071 ) compared many estimation methods for local galaxies 
and also showed the differences (the factors up to 2) among them. 

Except for above methods, the stellar mass could also be directly calculated using mid- 
infrared (MIR) luminosities in 3-5/zm, since emissions from the photosphere of old stars 
dominate the output of galaxies in this wavelength range. At the same time, the extinction 
and reddening in this range are clearly weak. But using these emissions to study the universe 
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was prevented because of the absorpti on of atmosphere until the onset of space infrared 
observations. Spitzer Space Telescope (IWerner et al.ll2004l ) is a very useful facility to help us 
to study the MIR to far-infrared (FIR) emission properties of galaxies. The two detectors in 
the shorter wavel ength bands (3.6/im and 4.5/xm) of Spitzer Infrared Array Camera (IRAC; 
Fazio et al.l 12004 ) could be treated as stellar mass tracers of nearby galaxies, although some 
disturbances may exist, such as the continuu m from hot dust or spectral features from 



the polycycl i c aromatic hydr ocarbons (PAHs; iLeger fc Puget I Il984c iPuget fc Leger I Il989 



Draine 


2003; 


Wu et al. 


2005) 



20031 : IWu et al.l 120051 ). Using about 150 local galaxies, iLi et al.l fl2007l ) showed that 
there truly existed tight correlations between the K-corrected 3.6/zm luminosities and the 
stellar mass based on the Bell et al.'s (2003) formula. 

In this paper, we want to derive the stellar mass formulae using IRAC 3.6 and 4.5 
jum luminosities based on the observation of Spitzer and Sloan Digital Sky Sur vey (SPSS; 
York et al. The Spitzer Wide-area Infrared Extragalactic Survey (SWIRE; lLonsdale et al 

20031 ). with a total field of ~49 deg 2 , is the largest extragalactic survey program among the 
six Spitzer cycle-1 Legacy Programs, and provides us a best opportunity to establish such 
correlations. 

The structure of the paper is as follows. We describe the construction of our sample and 
the referenced stellar mass in §2. The major results of the MIR stellar mass estimations are 
presented in §3. Discussion and summary are given in §4 and §5. Throughout this paper, 
we adopt a ACDM cosmology with Q m = 0.3, Q\ = 0.7 and H = 70kms~ 1 Mpc~ 1 . 



2. The Sample and Reference Stellar Mass 



Th e optical spectra l sample galaxies are selected from the S loan Digital Sky Survey 
(SDSS; lYork et al.l 120001 ) main galaxy sample ( Strauss et al.ll2002l ). with r-band Petrosian 
magnitudes lower than 17.77 mag. A total 15 deg 2 of SDSS sky matches the three northern 
Spitzer SWIRE fields. The derived data from SDSS DR7, including the stellar mass, was 
supplied by The Max-Planck-Institute for Astrophysics (MP A) and the Johns Hopkins Uni- 
versity (JHU) in their archives 0. Compared with former data releases, such as DR4, some 
galaxies were excluded in DR7 due to new restricts on the detected sources before publica- 
tion by MP A/ JHU . The stellar mass derived from fitting photometries with population 
synthesis models was treated as reference. The foreground extinction in all the five SDSS 
bands were corrected by subtracting the extinction values presented in the SDSS catalog. 
The IDL code by Blanton (version tag v4_l_4) was used to calculate the K-correction. The 



1 http://www. mpa-garching.mpg.de/SDSS/DR7 
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met hod and the SEP used in this code were described in details by iBlanton et al.l (120031 ) 
and IBlanton &: Roweid (120071 ) . 



The IRAC four bands (3.6, 4.5, 5.8 and 8.0 /im) images were mosaicked from the Ba- 
sic Calibrated Data (BCD; acquired from the Spitzer Scie nces Center) after the flat-field 



corre ctions, dark subtraction, linearity and flux calibrat i ons (IFazio et al. 



200J), with the final pixel scale of 0.6" flWu et al.ll2005l : ICao fc Wu 



2004 



2007 



Huang et al. 



Cao et al.ll200Rh . 



The MIPS 24/xm images were mosa i cked as the sam e way but with a pixel scale of 1.225" 



flWen et all 120071 : ICao fc Wu 1 12007 



Wu et al 



Skrutskie et al. 



20071). Based on the ca talogs of the Two 



19971 ; ICutri et al.l 120031 ) , we enhanced the 



Micron All Sky Survey (2MASS; 
accuracy of the astrometric calibration of 0.1" in all fi ve bands. The fin al IRAC and MIPS 
24 /im flux has calibration uncertainties less than 10% ( IRieke et al.ll2004f ). Then we matched 
these MIR sources with the above SDSS sample galaxies with a cross radius of 2". For the 
purpose of this work, the sources with both IRAC 3.6/zm and 4.5/zm detections were selected. 
The basic sample contains 1454 obje cts. For all these sources, a set of Spectral Energy Dis- 
tribution (SED) (IHuang et al.l 120071 ) was used to perform K-correction in IRAC 3.6/im and 
4.5/xm bands. 

Only part of the galaxies (total is 659) have been mapped by 2MASS. The Extended 
Source Catalog (XSC) was download from the archives of Infrared Processing and Analysis 
Center [IP AC] @ and matched with the above sample with a cross radius of 2". The K- 



correction for 2MASS NIR band fl ux was also based on the IDL code by IBlanton et al. 



(120031 ) and IBlanton fc Roweisl (120071 ) . The difference in K s ban d between Vega mag nitude 
which was used by 2MASS catalogs and AB magnitude is 1.84 (iFinlator et al.ll2000i ). And 
a factor of 5.12 (IFinlator et al.l |2000| ; iBinney fc Merrifield Ill998l )was used as the K s band 
solar absolute magnitude. 



2.1. Spectral Classification 



For all the 1454 galaxies, 68 with positive Ha equivalent width were classified as ab- 
sorption line galaxies, while the rest 1386 were classified as emission line galaxies. Futher- 
more, emission line galaxies with detected Ha, H/3, [OIII], [Nil] emission lines were se- 
lected, and optical spectral classifications were carried out on these 11 52 galaxies, adopting 



the traditional BPT diagnos tic diagram: [Nil] /Ha versus [OIII]/H/3 ([Baldwin et al. 



Veilleux &; Osterbrock ll987h . as sho wn in Figure [D The d ashed curve is from 



1981 



Kauffmann et al. 



(j2003bl ) and the dotted curve is from lKewley et al.l (120011 ). Objects located below the dashed 



2 http://irsa.ipac.caltech.edu/ 
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curve were classified as star-forming galax ies; those between these two lines w ere classified 
as composite (starburst + AGN) galaxies ( Kewley et al. 2006 ; Wu et al1ll998l ); while those 
above the dotted line were classified as narrow-line AGNs. Therefore, the final sample con- 
tains 1220 galaxies, including 561 star-forming galaxies, 292 composites, 299 narrow-line 
AGNs, and 68 absorption line galaxies, and 544 of this final sample have the 2MASS K s 
band flux. 



Wu et al.l ( 120071 ) presented an equation to compute (Iagn to represent AGN activities, 
defined as the distance of an AGN from its position ( x p ,y p ) in dex in the traditionally 



line-diagnostic diagram [OIII]/H/3 versus [Nil] /Ha to the lKewley et al.l (120011 ) 's curve along 
the parallel of the best linear fitting of 27 AGNs. This is quite similar to the star-forming 
distance defined by lKewley et al. (b006h in the diagnostic diagram [OIII]/H/3 versus [OI] /Ha 
or [SII]/Ha. Here, we also use this formula to quantify AGN activities (see the orange 
solid line in Figure [1] as an example) . In Figure CD, t he two dotted-d ashed lines above the 
dotted curve marked the boundaries that we re used bvlWu et al.l (120071 ). Only sources whose 
EW(H/3) were less than -5 were selected by IWu et al l (EoOTt ). But in this work, we did not 
use this limit to reduce those sources with smaller values of EW(H/3), such as LINERs. 
Therefore, many AGNs appeared outside the right boundary in Figure [TJ 



2.2. Define Referenced Stellar Mass 



The thermally pulsing asymptotic gi ant branch (TP-AGB ) could domina te the output 
of galaxies at some evolutionary phases (IMaraston et al.ll2006l ; iBruzual II2007I ). But lots of 
evolutionary tracks used in stellar population synthesis models did not consider the effect 
from TP-AGB. So there must exist bias for the derived stellar mass when comparing obser- 
vations wit h those models. For ex amp le, for young tem plate models (with ages less than 2 
Gyr), both IMaraston et al.l (120061 ) and IBruzual I (120071 ) have concluded that galactic stellar 
masses calculated with an improved treatment of TP-AGB stars are roughly 50 to 60 per 
cent lower. Hence, due to the lower effective surface temperature of TP-AGB, referenced 
stellar mass derived from optical observations could avoid the emission of TP-AGB. 

The referenced stellar mass used in this work was provided b y MPA/JHU in their 
archive %. They derived stellar masses from fitting photomet ries ( Salim et al. 2007 ) and 



spectral feat ures dGallazzi et al. 



from BC03 (IBruzual fc Chariot 



20051: iKauffmann et al.ll2003al ) with a large grid of models 



20031 ) spanning a large range in SFH. The stellar mass 



derived using SDSS photometries could avoid the aperture effect which should be considered 



3 http://www. mpa-garching.mpg.de/SDSS/ 
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when fitting spectral features, since only a part of region of the observed object was covered 
by SDSS spectral fiber with 3" diameter. For about eighteen thousand galaxies randomly 
selected from MPA SDSS DR4 Sample, the stellar masses derived from photometries and 
spectral indices fits were compared and shown in Figure 0(a). We found a good correlation 
between them. 

If the luminosity in some wavelength ranges could be obtained, we can estimate the 
stellar ma s s by a ssuming the mass-to- luminosity ratio in corresponding wavelength range. 
Bell et al.l ( 120031 ) presented the relationship between various optical colors and mass-to- 
luminosity ratios. In Figure E[b ), we compared t he referenced stellar mass with those com- 
puted based on t he equation bv lBell et al.l (|20031) . Here, to correct the 'diet Salpeter' IMF 



Bell et al 



f lBell et al 



2003) to the normal Salpeter (ISalpeter I Il955l ) IMF, we added a factor of 0.15 



20031 ) to the stellar mass calc ulated by using rest-frame g-r color and r band 



model magnitude. And a factor of 4.67 (IBell et al.l 120031 ) was treated as the r band solar 
absolute magn i tude (note there is a difference with a factor of 0.09 against the value used b y 
Blanton et al.l (120031 )). For the stellar mass derived from the equations of IBell et al.l ( 120031 ). 
there is a systematical overestimation with factor of about 0.2—0.6 dex, and such discrepancy 
is more obvious for a galaxy with lower stellar mass. This trend is more clear in Figure [3J 
The disagreement is more obvious for lower stellar mass and bluer ga laxies than higher stellar 



mass and redder ones, which is similar to the results supplied by of iKannappan fc Gawiser 



(120071 ) who found that large discrepancy existed between their reference stellar mass and the 
one derived from Bell et al.'s (2003) correlations, with factors in the range of 2 to 5. 



2.3. Sample Distribution 



The distributions of the SDSS B-band absolute magnitude (Mb), redshift, referenced 
stellar mass and r band concentration parameters of star-forming galaxies, composite galax- 
ies, AGNs and absorption line galaxies are shown in Figure SI Here, Mr were calculated 
from the SDSS g and r-band model magnitudes according to lsmith et al. ( 2002 ). Most of 
the galaxies, especially for composites, AGNs and absorption line galaxies, have Mb brighter 
than —18 mag, which us ually is regarded as the boundary to distinguish with dwarf galaxies 
( iThuan fc Martin 1119811 ). Few galaxies have redshift larger than 0.15. Herce, most of the 
galaxies in our sample are local ones. The star-forming galaxies have lower stellar mass than 
the others in Figure H(c). The morphology of galaxies was represented by r band concentra- 
tion parameters (R50/R90), which defined as ratio of two radii containing 50% and 90% of 
the Petrosian r band luminosity. The values of R50/R90 are smaller for early- type galaxies 
than those of late-type galaxies. According the distributions of four different spectral type 
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galaxies in Figure H^d), we found that the morphologies of the composite galaxies, AGNs 
and absorption line galaxies are similar, but those of star-forming galaxies are later type. 
Therefore, on the basis of above discussions, we find that low luminous dwarf early-type 
galaxies are scarce in this sample. 



3. Result: Stellar Mass Estimated by MIR luminosities 



A lot of factors could effect the stellar mass estimations. Therefore, good stellar mass 
tracers should be insensitive to these factors. Figure |5] show the correlations between mass-to- 
luminosity ratios of 2MASS K s band and two MIR bands luminosities and various factors. 
The K s luminosity was proved to be an ideal stellar mass tracer. Figure |5^al, a2, a3) present 
the effect of EW(Ha) on the mass-to- luminosity ratios. EWfHcQ is Ha equivalent width and 
could be regarded as the representative of SFH (IKennicuttl 1 1 9981 ) . The concentration param- 
eters (-R50/.R90) in Figure EJ^bl, b2, b3) represent the morphology of galaxies. The effect 
of color excess E(B-V) on the mass-to-luminosity ratios is illustrated in Figure E^cl, c2, 
c3). E(B - V) is for the intrinsic extinction and calculated from Balmer decrement F& a /F&p 
(ICalzettil 1200 ll ). Thus only the emission line galaxies are plotted here. In Figure [5](dl, 
d2, d3), the oxygen abun dance is used to represent the metallicity of star-f orming galaxies 
( Kauffmann et al.ll2003al ). The oxygen abundances of ISM were provided by lTremonti et al. 
(120041 ). which computed the abundances by fitting serial spectral features from SDSS obser- 
vational spectra to models. The effect of AGN act ivities is shown in Figure [5^el, e2, e3) 
by using cIagn to represent AGN activiti es as like Wu et al.l (120071 ) . Only AGNs located 
in between the two boundaries defined by IWu et al.l (120071 ) were adopted. We should note, 
except the concentration parameters, all above factors are derived from the spectral features 
in central regions covered by SDSS fibers with 3" diameter. For the correlations between 
M/L(3.6/im) and the five factors (EW(Ha), R50/R90, E(B-V), oxygen abundance, cIagn), 
the Spearman Rank-order correlation analysis coefficients are 0.74, 0.58, 0.15, 0.45 and 0.01, 
respectively; while for the correlations between M/L(4.5/xm) and the five factors (bottom five 
panels in Figure [5]), the Spearman Rank-order correlation analysis coefficients are 0.77, 0.58, 
0.10, 0.36 and 0.01, respectively. Hence, it is clear that the M/L(3.6/im) and M/L(4.5/im) 
are more sensitive to the EW(Ha), morphology and metallicity of galaxies. The value of the 
EW(Ha) represents SFH for galaxies without strong nuclear activities, while the morphology 
and metallicity are sensitive to the nowaday SFR other than the SFH. 

In order to derive the formulae to estimate stellar masses using these three band lumi- 
nosities, we present the correlations between them and stellar masses for all the galaxies in 
our sample firstly. The K s band, 3.6/im and 4.5/xm bands luminosities are plotted against 
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stellar masses in Figure [6j The stellar masses are the referenced one described in §2.2. These 
galaxies show good correlations between the three band luminosities and stellar masses. The 
best nonlinear and linear fits are shown as solid and dotted lines in Figure El and the fit- 
ting parameters are listed in Table 1. Here, the best nonlinear fits are obtained by using 
two-variable regression. Similar to Li07, we can derive star mass formulae based on the 
non-linear correlations: 



Log 10 M^f = (-1.60 ± 0.05) + (1.12 ± 0.02) x Log w vL u [K b ] 

Log l0 M^ ref = (-0.79 ± 0.03) + (1.19 ± 0.01) x Log 10 ^[3.6/im] 
Log w = (-0.25 ± 0.03) + (1.15 ± 0.01) x Log w uL u [A.5fim] 



(1) 

(2) 
(3) 



The correlations between two IRAC band luminosities and referenced stellar mass are as 
tight as that between the K s band luminosities and referenced stellar masses. The fitting 
residuals' standard deviations and the Spearman Rank-order correlation analysis coefficients 
of the correlations of 3.6/im and 4.5/xm luminosities with stellar masses are 0.11 and 0.12, 0.96 
and 0.95, respectively. The dashed line shown in Figure E^b) represents the non-linear fits for 



local luminous galaxies in the northern Spitzer SWIRE fields obtained by iLi et all (120071 ). 
Comparing the best non-linear fits in this work with that of ILi et all (120071 ). an obvious 
downward shift exists, especially for the non-luminous ones. For example, for galaxies with 
3.6/im luminosity of 1O 9 L , the stellar mass estimated fr om formula (2) would be about 
0.5 dex lower than the result from corresponding formula in lLi et all (120071 ). The shift could 
be due to fact that different stellar mass references (Figure [2]) are used in these two works. 

Above analysis ignored the effect by different SFH. In fact, there is such effect, at least 
for star-forming galaxies illustrated in Figure [5](al, a2, a3), which showed the correlations 
between the EW(Ha) and mass-to-luminosity ratios of K s band, 3.6/im and 4.5/im luminosi- 



ties. T he S FH could also be represent ed with the g-r color indicated by both lBell fc de Jong 



(boOlh and iKauffmann et all (l2003ah. which wa s also proved an effective method to esti- 
mate the stellar mass by iGallazzi fc Bell I (120091 ) . Figure [7] shows the correlation between 
g-r color and the EW(Ha) of our sample galaxies. There exists obvious anti-correlation 
between -EW(Ha) and g-r color. Therefore, we calibrate the correlations between mass-to- 
luminosity ratios and g-r color. Panels in Figure [H] show the correlations between various 
mass-to-luminosity ratios and g-r colors. Here, SDSS model magnitudes are used, and the 
foreground extinction and redshift have been corrected. 

These mass-to-luminosity ratios that are shown in Figure [H] are not constantly following 
the changing of color. From these four panels, we find out that redder galaxies tend to 
have larger mass-to-luminosity ratios. Using the two-variable regression, we obtain the best 
non-linear fits illustrated as solid lines in four panels of Figure [8] and the fitting parameters 
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are listed in Table 2. According to these correlations, we derive the formulae to estimate the 
stellar mass: 

Logio—rr^ = (-°- 81 ± °- 01 ) + ( L47 ± °- 01 ) x (9 - r ) ( 4 ) 



Log 



10 



Log 



10 



vL v {r) 

M*£ef 

vL v {K s ) 



Log 



in 



i/L„(3.6) 
uL v {A.5) 



-1.29 ± 0.05) + (1.42 ± 0.06) x (g - r) 
(0.23 ± 0.01) + (1.14 ± 0.01) x{g-r) 
(0.39 ± 0.01) + (1.17 ± 0.02) x(g-r) 



(5) 
(6) 
(7) 



The dotted-dashed line in Figure [5(a) is the respective line by iBell et al.l (120031 ) (see their 
Table. 7, Fig. 6 and Fig . 20). T he slope of our fit of M /L(r) vs. g-r color is much deeper 



than that of IBell et al.l (120031 ). iGallazzi &: Bell I (|2009[ ) also found similar discrepancy. One 
possible explanation is that different evolutionary population synthesis models have been 
used to derive standard stellar mass. The correlation between M/L(r) and g-r color is 
the tightest one, with the Spearman Rank-order correlation analysis coefficients of 0.97 
and the fitting residuals' standard deviations of 0.03. The correlations of M/L(3.6/zm) 
and M/L(4.5/xm) with g-r color are also tight, with fitting residuals' standard deviations 
of 0.05 and 0.06, and the Spearman Rank-order correlation analysis coefficients of 0.83 
and 0.80, respectively. The scatter in Figure E(a) is apparently much smaller than those 
in other panels, but the slope of the correlation between M/L(r) and g-r color for red 
galaxies is a little shallow er than the slope for blue ones. This variation is also suggested in 



Kauffmann et al.l (j2003al ) based on about 100,000 SDSS galaxies. Besides SFH, some other 



factors could possibly account for this slope variation, too. This variation can not be found 
in other panels in Figure El If the variations exist in Figure MJo, c, d), we can not clearly 
see them due to large scatters. The distribution of mass ratios of derived stellar mass and 
the reference one in Figure [9] show that the stellar masses derived by using Eq.6 and 7 are 
much tighter than the respective ones from Eq.2 and 3. 



4. Discussion 

4.1. Comparion with K s band 

K s band luminosity was often used to derive stellar masses of galaxies. From Figure El 
Figure [6] and Figure [5], we find that IRAC 3.6/im and 4.5/im luminosities are also good 
stellar mass tracers compared with K s band luminosity, because the radiation in all these 
bands was dominated by old stellar populations. Additionally, for the same sky area, the 
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effect of extinctions in IRAC 3. 6Atm and 4.5£un photometries is absolutely weaker than 



that in 2MASS K s band (e.g., iGao et al.l 120091 ) . There has been very few calibrations 



between stellar masses and 3/zm to 5/im emissions. One major reason is the limit of the 
weak transmissions of earth's atmosphere in this wavelength range. Now, the progress of 
space astronomy make us capable to observe the MIR sky directly. 



Bell et al.l (120031 ) indicated that the dispersion of blue galaxies was diffused in the 



M/L(K S ) vs. B-R plane (their Figure 20). And the fitting slope was shallower for red galax- 
ies compared to the blue ones. Using modified broad band Johnson photometry correlation, 



Kannappan fc Gawiser I (120071) also found sim i lar va riation of the slope. The transformation 



point of the slope of iKannappan fc Gawiser I (120071 ) was 1.2 in Johnson optical color B-R, 
which co rresponded to a SPSS g-r color of 0.55 based on the transformation methods sup- 
plied by iFukugita et al.l (119961 ). In their work, the slope of the correlation of M/L(K S ) vs. 
B-R was 0.5 for blue galaxies with B-R color l ess th an 1.2, while for the red galaxies with 
.B-the corresponding slope was 0.34. iBell et al.l ( 120031 ) found out that the reason for diffused 
dispersion and transformation of slope was the difference of metallicity of blue galaxies. For 
two galaxies with same in M/L(K S ), the one with lower metallicity must be bluer than the 
higher metallicity one optical colors. However we do not find such phenomena in other three 
panels( Figure [8]^b,c,d)), because of the lack of low metallicity galaxies in this sample, or the 
larger dispersion. Therefore, it must be more careful when estimating stellar masses of low 
metallicity galaxies using NIR or MIR luminosities. We will investigate their properties in 
more details in future. 



4.2. Effect of Star Formation 

Besides the EW(Hct), mass-to-luminosity ratios of two MIR luminosities are also sensi- 
tive to the concentration parameters which could be seen in Figure [5j Figure [10] shows the 
correlations between stellar mass ratios (the stellar mass derived from above formulae di- 
vided by the reference stellar mass) and morphology (R50/R90). The derived stellar masses 
in Figure ITuT a. b) base on nonlinear fits of 3.6/im and 4.5/im luminosities (with superscript 
'A') vs. reference stellar masses (Eq.2 and 3); while those in Figure [TuT c. d) base on nonlin- 
ear fits of M/L(3.6/im) and M/L(4.5/im) vs. g-r color (Eq.6 and 7) (with superscript 'B'). 
We find out that the obvious discrepancies of derived stellar masses of galaxies with different 
morphologies in Figure [TuTa. b) would overestimate the stellar masses of the late- type galax- 
ies compared to the early-type ones; while the ratios keep nearly constant in Figure [TuTc. 
d). SFH could account for the discrepancies in Figure fTUT a. b) when comparing the slopes 
of two upper panels and two bottom panels. Except for SFH, SFR may be another factor 
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resulting in overestimations of stellar masses for late-type galaxies. Galaxies with different 
morphologies have different SFRs . There are more new stars born in late- type galaxies than 
in early-type ones ( ILi et al.l 120071 ). SFR could relate with the values of intrinsic extinction 
and metallicity as shown in Figure Hence, it is necessary to check SFR on our stellar mass 
formulae. 

Panels in Figure [11] are similar to those in Figure IT0| just replacing r band concentration 
parameters by 24/im luminosities of galaxies. Same as Figure [Tnj stellar masses in two 
upper panels (Figure fTTTa. b)) are derived by nonlinear fits of 3.6/zm and 4.5/im luminosities 
(Eq.2 and 3), while those in bottom panels (Figure [TlT c. d)) based on the nonlinear fits of 
M/L(3.6/mi) and M/L(4.5/im) vs. g-r color (Eq.6 and 7 ). A template SED of a normal HII 
galaxy NGC 3351 (from SINGS; fcennicutt et al.l 12003^ was used to perform K-correction 
for the 24/im band for all the sample galaxies with 24/im detection. Figure [TT1 shows that Eq.2 
and 3 would overestimate the stellar mass for galaxies with higher 24/xm luminos ities. The 



24/im luminosity was proved to be a good SFR tracer by Spitzer observations (jWu et al. 



Calzetti et al. 


2005. 


2007; 


Zhu et al. 


2008) 



galaxies with lower star-forming activities. While Figure [IHc, d) are almost flat. Thus, the 
Eq.6 and 7 are not quite sensitive to SFR or morphology (Figure [TU]) . 

In Figure [TT], there are three outliers which are marked with red crosses. The one with 
number '3' is an AGN, which will be discussed in next subsection. Other two sources are star- 
forming galaxies with apparent features of interactions or mergers in SDSS images. The one 
with number '1' is a dwarf star-forming galaxies, with low metallicity of 7.93 and very large 
24/im-to-3.6/iHi flux ratio. The one with number '2' is a luminous star-forming galaxies with 
strong PAH emissions in IRAC 8/un band. Both of them have very blue g-r color. Actually, 
they are the only two galaxies with g-r color less than 0.1 in Figure [HJ Additionally, the 
stellar mass of '2' derived by using SDSS photometries and spectral indices are apparently 
different. Therefore, it is hard to define which one is the best estimator. The stellar mass 
estimated using photometries is about 1 dex less than the one estimated by spectral indices. 
For '1' and '3', there are no stellar masses values derived by spectral indices. 



Apart from emission directly from the photosphere of old red stars, t 



re continuum from 



very small dust hea ted by hot stars, even some broad emissio n features (Gillett et al. 19731; 



200< 



3; ICao fc 



Wu 



Will ner et al.l 119771 ) which were indicated to be from PAHs (I Wen et al. 
20071 ). could increase the MIR luminosity of a galaxy. The dust emission dominates radiation 
of star-for ming galaxies in th e wavelength range beyond 5/im, but could also enhance the 
4.5/im flux tlHelou et al.l 120041 ) . To evaluate such effect, the 4.5/im-to-3.6/im flux ratio is 
plotted against the 24/im luminosity in Figure [12] The three outliers in Figure [11] were 
marked with red crosses. The blue boxes are the AGNs. There are some AGNs with very 
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high 4.5/im-to-3.6/im luminosity ratios marked with red boxes. Except these outliers, we 
find a correlation between 4.5/im-to-3.6/im luminosity ratio and star formation for most of 
sources. Because the 24/im luminosity is the tracer of SFR, this correlation may be due to 
the dust emissions in the 4.5/im band. With the assumption that the entire IR AC 3.6/im 
band luminosity is from the stellar emission, a factor of 0.596 (IHelou et al.l l2004j ) was used 
to scal e the stellar continuu m of 3.6/im to that of 4.5/xm based on Starburst99 synthesis 
model (ILeitherer et al.l Il999f ) with the solar metallicity and a Salpeter initial mass function 
(IMF) between 0.1 and 120 This factor is similar to that derived from the observed 

SED of early type galaxies by IWu et al.l (120051 ) . The factor of 0.596 is plotted with dashed 
line in Figure [121 If this factor is available for all star-forming galaxies, we could estimate 
the dust emi ssion fraction ba sed on the correlations between 24/im and total-infrared (TIR) 
luminosities (lZhu et al.l 120081 ). e.g., for luminous infrared galaxies with TIR luminosity of 
1O 11 L , about 30% flux at rest-frame 4.5/xm is from dust emission. Therefore, compared 
to the 4.5/im emission, IRAC 3.6/im emission is a better stellar mass tracer. This is also 
supported by the much tighter correlation between 3.6/im and stellar masses in Table 1 and 
2. 



4.3. Effect of AGNs 



Due to the powerful UV emissions from AGNs, the dust surroundin g them could be 



heated by central monste r s and then re-radiate in IR range (see, e.g, iBell et al.l 12005 
Perez- Gonzalez et al.l lioiiS felbaz et~al1l2007l : iDaddi et alJl2007h . Higher MIR-to-H a/MIR- 
to-UV luminosity ra t ios were found b y previous works based on Spitzer observations (IWu et al 
20071 ; iLi et al.l 120071 : IZhu et al.l 120081 ) . So we want to know whether the energy from AGNs 
could strengthen the 3.6/im and 4.5/im emissions. Figure [£3] shows the correlations be- 
tween stellar mass rati os and AGN act ivities. Here, j ust the AGNs located in side the two 
boundaries defined by IWu et al.l (120071 ) are adopted. iKauffmann et al.l (l2003bl ) found that 
the contribution from AGN light to the optical continuum was much weaker compared to 
from starburst. Thus, ignoring the effect of AGN emission on the estimation of referenced 
stellar mass would be reasonable. Red boxes represent the re spective AGNs in Figure [12] 
with 4.5/im-to-3.6/iHi flux ratio large than 0.9. According to IWen et al.l (120071 ). we know 
that AGNs with very high 4.5/im-to-3.6/im ratios may be QSOs. In Figure [131 except those 
scatters, the stellar mass derived from MIR luminosities keep constant with increasing cIagn, 
while Figure [5](e2, e3) show that AGN activities can not disturb mass-to- luminosity ratio 
for AGN-hosting galaxies, even at 4.5/im(Figure fT4|) . Therefore, there seems no confirmed 
non-thermal radiations or dust emissions heated by AGNs in IRAC 3.6/im and 4.5/im de- 
tections. Because all the sources here are type-II AGNs, the powerful intrinsic extinction of 
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dust torus of narrow line AGNs could be another explanation for non-detection of AGNs in 
3.6/im and 4.5/zm in our sample. The absence of AGN emissions at 3.6/zm and 4.5/xm also 
indicate that the stellar mass formulae based on 3.6/xm and 4.5/zm luminosities are applicate 
to most AGN-hosting galaxies. 

Figure 0] present the distributions of absolute B-band magnitudes, reference stellar 
masses and r band concentration parameters. The distributions in Figure |4] of composite 
galaxies and AGNs are strikingly different from those of star-forming galaxies, but similar to 
absorption line galaxi es. These results indicate that AGNs in o ur sample tend to be hosted 
in early-type galaxies (jKauffmann et al.ll2003bl ; IZhu et al.ll2008l ). Based on these results and 
the discussions about the effect of star formation in §4.2, we find star formation rates in our 
AGN-hosting galaxies are lower. 



5. Summary 

Based on the sample cross-identified from Spitzer SWIRE field and SDSS spectrographic 
survey, we derived the formulae to calculate the stellar mass using non-linear correlations 
wth IRAC 3.6/i and 4.5/xm luminosities. The mass-to-luminosity ratio of these two MIR 
luminosities are sensitive to the EW(Ha), morphology and metallicity of galaxies, especially 
for the former one, while the value of the EW(Ha) represents the SFH for normal galaxies. 
Using the g-r color to represent SFH, we re-calibrate the stellar mass formulae for various 
galaxies. We find these formulae are better than those which ignore the effect of SFH. We 
could not conclude the applicability to low metallicity galaxies due to a lack of such object 
in our sample, thus one must be really careful when adopting our formulae to estimate 
stellar masses of those galaxies. Additionally, we found that the dust emission heated by 
hottest young stars could enhance the IRAC 4.5 /im luminosity, especially for luminous 
infrared galaxies. So the formulae we derived are not applicate to these galaxies. And these 
formulae are also not applicate to the stronger AGNs, such as QSOs, although almost all the 
AGN-hosting galaxies in our sample present similar properties like absorption line galaxies 
(early-type), with low SFRs. 
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Logio [Nll]/Ha 



Fig. 1. — The BPT di agnos tic diagram: [ Nil] /Ha vs [0III]/H/3. The criteria from 



Kauffmann et all ( 12003bl ) and iKewley et all ( 120011 ) are illustrated as dotted and dashed 



curves, respectively. The objects below the dotted curve are defined as star-forming galax- 
ies. The triangles between the above two curves are those classified as composite galaxies, 
while the boxe s above the dash ed curve denote AGNs. We use the distance cIagn, which 
was defined by lWu et all (120071 ). to ch aracterize AGN a ctivity. The two dotted-dashed lines 
mark the boundaries that were used by IWu et al.l (120071 ) . An example was show in this figure 
with orange solid line. 
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Fig. 2. — Comparing the referenced stellar mass, which derived by photometries fitting, 
with those from spectral features fitting supplied by MPA/JHU (Panel (a)) for about 
eighty thousand galaxies randomly selected from SDSS DR4. In Panel (b), the comparing 
between the reference and the stellar mass computed based on the rest-frame g-r color, r 
band model magnitude and the equation bv lBell et al.l (20031) is presented. Please note that 
a factor of 0.15 was adding t o the equation o f iBell et al.l (120031 ) to correct the 'diet Salpeter' 
IMF to the normal Salpeter jSalpeter Wl95m IMF. 



-20 - 



Fig. 3. — The correlations between the stellar mass ratio and the referenced stellar mass 
(a), g-r color for about eighty thousand galaxies randomly selected from SDSS DR4. The 
stellar mass ratio represent the discrepancies between the referenced stellar mass and the 
stellar mass estimated based on the rest-frame g-r color, r band model magnitude and the 
equation by iBell et al. is pr esented. Same as in Figure [2j a factor of 0.15 was adding 

to the equation of lBell et al.l ( 120031 ) to correct the 'diet Salpeter' IMF to the normal Salpeter 
jSalpeter Hl95fi( l IMF. 
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Fig. 4. — The distributions of (a) absolute B-band magnitude; (b) redshift; (c) the referenced 
stellar mass; (d) the concentration parameters (R50/R90) for star forming galaxies (black 
solid lines), composites (red dashed lines), AGNs (blue dotted lines) and absorption line 
galaxies (pink dashed-dotted lines). 
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Fig. 5. — Comparing the mass-to-luminosity ratios with various factors, include Ha equiv- 
alent width (al, a2, a3), R50/R90 (bl, b2, b3), intrinsic reddening (cl, c2, c3), oxygen 
metallicity (dl, d2, d3) and the activity of AGN (el, e2, e3). The 2MASS K s band lumi- 
nosity is used in upper five panels. The 3.6/zm luminosity is used in middle five panels. The 
4.5/xm luminosity is used in the bottom five panels. 
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Fig. 6. — Correlations between referenced stellar mass and 2MASS K s band, two IRAC 
bands luminosities for all the galaxies. The best nonlinear and linear fits are illustrated as 
solid and dotted lines. The red dashed lines show in Panels ( b) represent the nonlinear fits 
for more than one hundred luminous SWIRE-field galaxies by iLi et al.l (120071 ). 
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Fig. 7. — Correlations between EW(Ha) and g-r color for all the galaxies. 
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Fig. 8. — Correlations between g-r color and various mass-to-luminosity ratios for all the 
galaxies. Due to the relatively lower sensitivity in 2MASS K s band observation, only 544 
galaxies were detected and presented in Panel (b). The solid lines r epresent the best nonlinear 
fits. The dotted-dashed line in Panel (a) is the respective line by iBell et al.l (120031 ). 
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Fig. 9. — The distributions of ratios of the derived stellar masses divided by the referenced 
ones. The derived stellar masses were computed based on Eq.2 (black solid lines in Panel 
(a)), Eq.6 (red dashed lines in Panel (a)), Eq.3 (black solid lines in Panel (b)), Eq.7 (red 
dashed lines in Panel (b)) 
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Fig. 10. — Correlations between concentration parameters and the stellar mass ratios, which 
are the discrepancies between the stellar mass derived from our calibrated formulae and the 
reference. The derived stellar mass by using Eq.2, 3, 6, 7 is shown in Panel (a), (b), (c), (d), 
respectively. 
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Fig. 11. — Similar as Figure 10, but the concentration parameters were instead of 24/zm 
luminosities. The red crosses represent the three scatters, and they were signed out with 
numbers. 
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Fig. 12. — Correlations between 24/im luminosity and the flux ratio between 3.6/zm and 
4.5/xm for all the galaxies. The blue boxes represent AGNs. The red boxes represent AGNs 
with 3.6 / um-to-4.5/im ratio larger than 0.9. The three outliers in Figure 11 were sighed out 
with red crosses. The dash ed line represent t he factor of 0.596 which is the ratio between 
4.5/xm and 3.6/im flux by (IHelou et al.l 120041 ) with the assumption that the entire IRAC 
3.6/iHi and 4.5/zm band luminosities are from the stellar emission. 
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Fig. 13. — Correlations between AGN activity and the stellar mass ratios, which are the 
discrepancies between the stellar mass derived from our calibrated f ormulae and the refer- 
ence. Here, the distance cIagn, which was defined by lWu et al.l ( 120071 ). to characterize AGN 
activity. The derived stellar mass by using Eq.2, 3, 6, 7 is shown in Panel (a), (b), (c), 
(d), respectively. The red boxes represent the respective AGNs in Figure [12] with higher 
4.5/im-to-3.6/zm flux ratio (large than 0.9). 
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Fig. 14. — Correlations between the distance cIagn, which represent AGN activity, and the 
flux ratio between 3.6/im and 4.5/xm for AGNs. The red boxes represent the sources with 
4.5/zm-to-3.6/im flux ratio large than 0.9. The dashed line represent the factor of 0.596 such 
as plotted in Figure [121 
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Table 1. Correlations coefficients between stellar mass and K s band, 3.6/im and 4.5//m 

luminosities 





X 


a 


6 


s 


r 


c 


N 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


M* 




- 1.60±0.05 


1.12±0.02 


0.10 


0.94 


0.28±0.11 


544 


M* 


^L„[3.6/ira] 


- 0.79±0.03 


1.19±0.01 


0.11 


0.96 


0.99±0.14 


1220 


M* 


uL u [4:.5^m] 


- 0.25±0.03 


1.15±0.01 


0.12 


0.95 


1.17±0.14 


1220 



Note. - Col.(l): the referenced stellar mass in solar unit; Col. (2): names 
of MIR luminosities; Col.(3)-(4): the coefficients a and b of the nonlinear fit: 
log 10 (y) =a+blog 10 (x); Col. (5): the standard deviation s of the fitting residuals; 
Col. (6): the coefficient r of the Spearman Rank-order correlation analysis; Col. (7): 
the coefficient c of the linear fit: log 10 (y) =c+ log 10 (a;); Col. (8): the number of 
sample galaxies used for the fitting procedures. 
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Table 2. Correlations of M*/L and colors 



y 


X 


a 


b 


s 


r 


N 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


M*/L v {r) 


g-r 


-0.81±0.01 


1.47±0.01 


0.03 


0.97 


1215 


M,/L U {K S ) 


g-r 


-1.29±0.05 


1.42±0.06 


0.08 


0.54 


410 


M^/vL l/ [3.6fj,m] 


g-r 


0.23±0.01 


1.14±0.01 


0.05 


0.83 


1210 




g-r 


0.39±0.01 


1.17±0.02 


0.06 


0.80 


1208 



Note. - Col.(l): various mass-to-luminosity ratios in solar unit; 
Col. (2): g-r color; Col. (3)- (4): the coefficients a and b of the nonlinear 
fit: log 10 (y) =a+blog 10 (x); Col. (5): the standard deviation s of the fitting 
residuals; Col. (6): the coefficient r of the Spearman Rank-order correlation 
analysis; Col. (7): the number of galaxies used for the fitting procedures. 



